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Thompson scattering of cosmic microwave background (CMB) photons off of free electrons during 
the reionization epoch induces a correlation between the distribution of galaxies and the polarization 
pattern of the CMB, the magnitude of which is proportional to the quadrupole moment of radiation 
at the time of scattering. Since the quadrupole moment generated by gravitational waves (GWs) 
' gives rise to a different polarization pattern than that produced by scalar modes, one can put 

interesting constraints on the strength of GWs on large scales by cross-correlating the small scale 
galaxy distribution and CMB polarization. We use this method together with Fisher analysis to 
£\j ■ predict how well future surveys can measure the tensor-to-scalar ratio r. We find that with a future 

CMB experiment with detector noise Ap — 2 ^K-arcmin and a beam width Afwhm = 2' and a 
future galaxy survey with limiting magnitude / < 25.6 one can measure the tensor-to-scalar ratio 
with an error oy ~ 0.09. To measure r w 0.01, however, one needs Ap ~ 0.5 /iK-radian and 
Afwhm — 1'. We also investigate a few systematic effects, none of which turn out to add any biases 
to our estimators, but they increase the error bars by adding to the cosmic variance. The incomplete 
sky coverage has the most dramatic effect on our constraints on r for large sky cuts, with a reduction 
in signal-to-noise smaller than one would expect from the naive estimate ijs) oc /sky Specifically, 
. we find a degradation factor of /d og = 0.32 ± 0.01 for a sky cut of |6| > 10° (/sky = 0.83) and 

/dcg = 0.056 ± 0.004 for a sky cut of \b\ > 20° (/sky = 0.66). Nonetheless, given that our method 
has different systematics than the more conventional method of observing the large scale B modes 
. directly, it may be used as an important check in the case of a detection. 
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I. INTRODUCTION 

The possible detection of a primordial B mode signal in the polarization of the Cosmic Microwave Background 
(CMB) will certainly be of great importance for cosmologists and high energy physicists alike as it will provide us 
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, with information about physics at high energies that will not be accessible through terrestrial experiments in the 
foreseeable future (see e.g. [l|). 

The leading mechanism for setting up the initial condition of the Universe is inflation 043) i m which the quantum 
fluctuations of a scalar field, the inflaton, inside the horizon gets stretched out of the horizon during an almost 
| exponentially expanding phase of the universe, generating the primordial seeds of structure in the Universe 043- 
T — I . The predictions of this theory have so far passed all the existing observational tests (see e.g. H). However, little is 
known about the properties of the inflaton field(s) and its potential. A tremendous amount of insight will be gained 
. £h . if the primordial B pattern of the polarization of the CMB can be measured. This pattern can only be generate d by 
primordial tensor perturbations, that is the gravitational waves (GWs) produced during the inflationary era 0, llfjj |. 
The spectrum of GWs is commonly expressed as A^(fc) = Ah(k*)(j-) nT , where is the power spectrum of the 
traceless-symmetric part of the metric, hij, per logarithmic interval in wavenumber, fc, and fe* is an arbitrary pivot 
scale. Furthermore, the strength of the tensor modes is commonly quoted in terms of the "tensor-to-scalar ratio" 
r = ^(fc*) i where Ar is the power per logarithmic interval in wavenumber in the curvature perturbations 1Z. The 
magnitude of r depends on the Hubble scale during inflation, which is in turn a function of the inflaton potential 
energy during inflation. Of special importance is the value of r ~ 0.01 since this value corresponds to a GUT scale 
energy during inflation so a detection of the signal at this value would strongly suggest a relationship between inflation 
and GUT scale physics. 

However, the signal in the B mode is very small and is peaked on large angular scales where the galactic foregrounds 
are most difficult to remove. In addition, the lensing of the CMB by large scale structure contaminates the signal by 
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transferring power from the E mode to the B mode [37jj . This contaminant can however be removed sufficiently by 
"delensing" techniques [H10] to reach the critical accuracy of r ~ 0.01 fl3l. IbH]. 

Having different methods of measuring the GW signal with different systematic errors is very important to make 
sure that the detected signal is primordial, not an artifact of instrumental or foreground contaminants such as lensing. 
The aim of this paper is to propose a new method to measure (or put a bound on) the strength of GWs by using 
the correlation between the galaxy distribution and the CMB polarization fluctuations. That there should be a cross- 
correlation is clear: Thompson scattering of CMB photons off of free electrons during the reionization era introduces 
anisotropies to the polarization pattern of the CMB [15| . The amplitude of these anisotropies is proportional to the 
number density of free electrons, which itself is correlated with the distribution of galaxies, hence the correlation 
between the galaxy distribution and the CMB polarization. Furthermore, these anisotropies are also proportional to 
the quadrupole moments of the CMB radiation at the time of scattering. The quadrupole moments that give rise 
to the B polarization pattern can only be generated by the tensor perturbation of the metric so by measuring the 
amplitude of the cross-correlation between galaxies and CMB polarization one get an estimate on the amplitude of 
GW signal. 

The outline of the paper is as follows. In section Hi Al we find an analytical formula for the cross-correlation between 
the CMB polarization patterns and the galaxy distribution. We will then use this result in section IIIBI to find a 
quadratic estimator to construct the electric and magnetic type polarization moments, E\ m and B\ m , generated by a 
smooth electron density field at a given redshift bin i from the observed galaxy distribution and the CMB polarization 
fluctuations. The details of this calculation can be found in Appendix [XJ The noise and signal covariance matrices for 
detecting these average polarizations are presented in section III CI and the derivation is presented in Appendix [B] In 
section HTT1 we use the Fisher formalism to forecast the power of this method to constrain the tensor-to-scalar ratio r 
for futuristic CMB and galaxy surveys. Several possible systematic effects are considered in section ITVl We conclude 
in section El 

WMAP5 parameters are assumed throughout the paper. 

II. FORMALISM 

A. Polarization from inhomogeneous reionization 

We define spin ±2 variables ± 2 P(n) on the sphere in terms of the Stokes parameters q and u as [To| 

± 2P(n) = {q±iu](n). (1) 
These can be found from the integrals of the temperature quadrupole along the line of sight, see e.g. (l7j . 
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P(n) = l^L J dDg(Dn) £ Q™(DA) ±2 Y 2>m (n), (2) 



m=-2 



where D is the conformal distance in units of the Hubble distance today, and g is the visibility function. 

If we only consider the polarization produced at a given redshift slice "i," and further separate the visibility function 
at that slice into a smooth part and a part coming from the fluctuations in the electron number density we will have 



g\Dn) = g\D) (^1 + j = ^(£)(i + AJ(Dn)) = g\D) + -i±—L j . ( 3 ) 

Here, A& and A g are the baryon and galaxy overdensities, respectively, and b is the galaxy bias with respect to the 
baryons. Similarly, the polarization generated at that redshift slice can be separated into a smooth part and overlying 
fluctuations. Then Eqs. © and (J3]) give 

±P l (n) = ± P J (n) + ( 5 ± P I (fi), (4) 

where 

l 2 

± P*(n) = l^L J dDg\D) £ Q^(Dh) ±2 F 2 , ro (A) (5) 



and 



S ± P\n)^ ±J r(n)A*(n)l (6) 
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Here A 1 is the projected overdensity of galaxies at that redshift bin 

A' (A) = f * r iSh dDAl(Dn). (7) 

To find the E and B polarization modes coming from that redshift slice we expand all of the functions in terms of 
their multipole moments 



We also use the identities 



1 1 II 

± P\n) = ±PL(±2Yi, m (A)), 

I III 

6 ± P\n) = J2 S ± P L(±2Y l , m (&)). (8) 



Elm — ~{+Pl m + -Plm), 

B lm = -(+Plm - -Plm)- (9) 
ll 



Then after a straightforward calculation we find 



SE L - EE( _1 ) m (-m mi m 2 ) F ^ u ^ ai ^ lE h^i ~ lhhi B i imi )— 



limx l 2 m 2 

where a, 7 and F are defined as 

_ ( Wl+fe+A 



SB L = E E(-!) m ( - m m\ m 2 ) WTfeW^U + OhW^U)^. (10) 



(11) 
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(13) 





1 




2< 




1 


TUfaJ = 






J 



(2Zi + 1)(2Z 2 + 1)(2Z + 1) ( I h I 



1-2 



4?r 12-20 



Estimators 



As we can see from the results of the previous section the power in SEi m /SBi m modes at a given I comes from 
a quadratic sum of a range of multipolcs in Ei imi / E>i imi and A/ 2 „ l2 , where I, l\ and I2 must satisfy the triangle 
inequalities. Since we are looking for large scale gravitational wave modes, i.e small li, and because the power in the 
galaxy distribution is large at small scales (large ^2), we realize that the triangles we are dealing with in I space are 
elongated with two long sides of length I and I2 and one small side of length l\. In other words, to find an estimate of 
the power of the CMB polarization anisotropics on large angular scales we can use the small scale power of the CMB 
polarization and the galaxy distribution. 

We show in Appendix [X] that an unbiased, minimum variance, quadratic estimator for E l j B % can be found to be 

eIm - ^fEEH) 1 ^; t 2 -m)^( l )K^<--^w<-) a ^^ (1 4 ) 



Zimi I21712 



hmi l 2 m 2 



bi m = ai b E E (-!) M ( 4 j? 2 - M ) U^iCu + OhhA) ( 15 ) 
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where 



9ff{L) = Ull l « . , (16) 

I1I2L I2 



M hl 2 L 
Kl 2 L 



(lau^Cf^ + l^ifC^), (17) 

(i7LMj 2 ^r obs +i^M 2 i 2 < s ' obs ), (is) 



flhh - F Lhh h 1 (19) 



A\ x = (2L+1) 



^2 fLhl 2 9hf 2 ( L ) 



(20) 



U1Z2 

Here X can be either E or B and £f^ s , BP^s and A^ s are the observed quantities which can be decomposed as 



J lm ' Im allKX t-^lm 

\robs 1 \^noisc yrccom , yreion 1 x^lcns 1 \^p.s. . x^gal , \^noisc 

A im — A im + A im — A Im + A im + A /m + A /m + A /m + A /m ) 

. i , Ai.no 



Ai,obs \i - a 2, noise foi "\ 



where the superscripts are for recombination, reionization, lensing, polarized point sources, polarized galactic fore- 
grounds and noise (detector noise in the case of polarization and Poisson noise for galaxy overdensities) . A quantity 
without any superscript description represents all the sources that contribute to it except the noise. The polarization 
from the reionization epoch can be separated into the contribution from different rcdshifts, 

XZ on = J2(^L + ^i m ), (22) 

i 

as are calculated above. The power spectra in the above equations are defined as 

(A&Af, m ,) = Cf gi 5 w 8 mml , 
(A ^bs* A i,ob S) = c? 3i > ohB 8 w 8 mm > = {Cf 9 ' +Nf a ')5 w 5 mm ,, 

( A im X Vm>) = C l "U'Omm' = [yi + )0U'd mm >. (23) 

The noise power spectrum for the galaxy-galaxy correlation is simply the Poisson noise, 

Nf a ' = \, (24) 

where n l g is the mean projected number density of galaxies at redshift bin "i" in units of sr _1 . 

The instrumental noise power spectrum for either E or B mode polarization can be written as [l8| : 

A f ^ 2 



V TcMB J 

where Ap is the detector noise in units of /xK-radian, X"cmb = 2.725 x 10 6 /iK and #fwhm is the width of the beam 
in units of radians. 

C. Signal and Noise covariance matrices for E % LM and B l LM 

The covariance matrices of our estimators for the average polarization generated at a given redshift bin can be 
written as 

{XI\Mm') = (Cf + Nf* 3 ) 5 W 6 MM ,, (26) 
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where C x x 3 and N x x 3 are the signal and the noise covariance matrices, respectively. The noise covariance can 
be calculated to be 

= Y.tii&* l2L C^. (27) 

Since we are looking at the small scale galaxy distribution it is justified to ignore the correlation between the galaxy 
distributions at different redshifts, i.e. Cf 9 — Cf 9 Sij. The formula for the noise covariance matrix then simplifies 
to 

A / V V ,1',/%. (28) 

Parity consideration or a direct calculation shows 

Nf S3 = 0. (29) 
The signal covariance matrix, C x x3 , for E and B type polarization can be written as 

Cf E3 = Cgf+Cgf, (30) 

Cf 53 = (tP • (31) 

where the subscripts S and T mean scalar and tensor, respectively. There is no contribution to the B polarization 
from scalar perturbations. It is worth mentioning that Cf g 3 actually acts as a source of noise for the purpose of 
detecting G Ws even though it is the main contributer to E\ m . 

We show in appendix [B] that the tensor terms in the above equations can be calculated as 

Cft? =^jdk k 2 S l { x {k) Sf' x (k)P h (k), (32) 
where Sf' X (k) is given in Eq. (|B12|) . We follow reference in writing the power spectrum of GWs as 

A 2 ^) = = Al(k ) ( A "I" (33) 



2vr 2 nx \ k 

for fco = 0.002 Mpc^ 1 . For the purpose of this paper we take the tensor spectral tilt to be zero, n t = 0. Furthermore, 
the amplitude of the GW spectrum can be parametrized in terms of the tensor-to-scaler ratio, r, as 



A 2 ,(fc ) 
A^(fc )- 



(34) 



where An(k) is the curvature perturbation spectrum. We fix A^(fco) = 2.41 x 10~ 9 . The tensor-to-scalar ratio, r, 
is then the only parameter in this paper that we try to constrain and the rest of the parameters are fixed to their 
fiducial values. 

To present the results of this section in a coherent way we show in Fig. [T]the steps one needs to take to estimate 
cosmological parameters using our method. 

III. SENSITIVITY TO GRAVITATIONAL WAVES 

In the absence of appropriate data we use the Fisher information method to forecast the capacity of a few futuristic 
surveys to detect GWs using our method. 

So far we have found expressions for the covariance matrices of the signal and noise. From them we define the total 
covariance matrix as 

CiLM-i'L'M' = (^l)u'SlL'SmM' , 

(\ L ) u ,=CE i3i ' +NE ij§i '. (35) 
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Measure the galaxy distribution 
at multiple redshift bins "\": Ai 



Use Eq. 15 to find an estimate for 
the polarization pattern generated 
by the smooth electron field at any 
redshift bi 



Do a likelihood analysis to determine the cosmological 
parameters. The likelihood from our method is 
-2 In L = In det C + x C 1 x' 
where the data vector x consists of B l LM , calculated in 
the previous step, and the covariance matrix C iLM . t i L i M r 
can be computed using Eqs. 26, 28 and 32 



FIG. 1: An schematic plot showing how our method can be used as a parameters estimation tool. 



Then assuming that both signal and noise are Gaussian random variables, the Fisher matrix can be written as 

Fab = ^Tr [C, B C- 1 C, 6 C- 1 ] , (36) 

where the derivatives a and b are with respect to the parameters we wish to constrain. The inverse of C can be 
written as 

(C )iLM;i'L'M' = (A L )**' $MM' ■ (37) 

Using this in Eq. (|36[) and performing the sums over the L and M indices we find 

O T _l 1 

Fab = £ ^T-Tr [(A L ), a (AZ 1 )(A L ), b (A Z 1 )] . (38) 

L 



The 1-sigma error bars on the parameter a marginalized over the other unknown parameters is y(F 1 ) aa , but 
assuming perfect knowledge of the other parameters this error is 1/ yj (F aa ). In our analysis we find constraints on 
the tensor-to-scalar ratio, r, by fixing all other cosmological parameters to their fiducial values so our Fisher matrix 
has only one entry and oy = \j\/F rr . 

To calculate the noise covariance matrix from Eqs. [TBTfSOl we need the polarization power spectra, C xx , the 

instrumental noise spectrum, N xx , the power spectra of the galaxy distribution at different redshifts, C? 9 , their 
noise spectra 9 and the bias b l . We calculate them as follows: 

We use the publicly available code CAMB 38] to calculate C xx including lensing. We use Eq. [55] to calculate 
N xx for several choices of detector noise A p and beam width 6*fwhm- These power spectra are shown in the left 
panel of figure [5] The noise power spectrum corresponds to a futuristic experiment with A p = 2 /iK-arcmin and 
^fwhm = 2' (for comparison, the detector noise for the ACTPol Deep survey is Ap = 4/xK-arcmin at v = 150 

GHz and that of SPTpol survey [Io| is A p = \J NEQ = 9 /iK-arcmin at 90 GHz. These are achieved over a 

small area but the technology is advancing rapidly and wide sky coverage could be feasible in the near future). The 
tensor-to-scalar ratio is put to its fiducial value of zero here, so the B mode is exclusively from the weak lensing of 
the primordial E mode. 

The noise in the galaxy power spectrum is Nf 9 — l/n l g , where n l g is the mean number of observed galaxies per 
steradian in the i th redshift bin. To find it we assume a survey similar to that of LSST with limiting magnitude of 
I < 25.6 together with a northern counterpart to have a full sky coverage. Our slices in redshift are in the range 
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FIG. 2: [Left] The polarization power spectra for EE (solid red) and BB (dashed green) modes from recombination + reionization 
+ lensing as computed by CAMB. The gravitational waves are not included so all the B modes are coming from the lensing of 
the primordial E modes. The noise power spectrum for a future experiment with A p = 5/iK-arcmin and #fwhm = 2 arcmin 
is also shown (dotted blue). [Right] Projected galaxy-galaxy signal and noise power spectra for a future survey with limiting 
magnitude of I < 25.6, shown for the first and last redshift bins centered at z — 0.5 and z = 4.1, respectively, both with a 
width of Az = 0.2. 



0.4 < z < 4.2 with a width of Az — 0.2. We show the noise power spectrum for the first and last redshift bins in 
the right panel of Fig. [5] together with the signals at the corresponding redshifts. The signals are calculated from the 
dark matter power spectrum by using the bias for the star forming galaxies 

b{z) = 0.9 + 0.4 z, (39) 

which is a fit to the results of Ref. [2l|. We assume baryons trace the dark matter distribution on large scales, so that 
the same equation for the bias can be used in Eq. 1191 

Given the above specifications we can calculate the one sigma error bars on r, oy = \/l/F rr . We show in Fig. 1 

the constraints one gets by constructing either the average magnetic type polarization at each redshift, B l LM: (red 

circles) or constructing the electric type polarization, E l LM , (green crosses) as a function of L. As expected, the B 
modes are more sensitive to the GWs since here the signal is not swamped by the scalar modes (see Eq. [50)1 . Also, 
since the signal drops rapidly with increasing L (see right panel of Fig. [5]) most of the information about the GWs is 
contained in its lowest moment, i.e. the quadrupole. 

From Fig. [5] it can be seen that the signal goes below the noise at different values of I for the E or B polarization 
modes and similarly for the galaxy power spectra at different redshift bins. An interesting question then arises: at 
what angular scales is most of the signal in determining Blm coming from? To find out, we have plotted in Fig. 0] 
the 1-sigma error bars on r as a function of the maximum I2 in the sum in Eq. [20] (the maximum l\ will automatically 
be determined from the triangle inequality and the value of L). We can see that oy decreases like a power law with 
increasing Z!f ax until ^ nax ~ 4000 where it starts to level off. 

To explore the dependence of oy on the specifications of a future CMB and galaxy survey, we show in Fig. [5] a 
contour plot of ay as a function of detector noise and beam width for two different redshift surveys. The contours are 
for 0> = 0.1, 0.05, 0.01 and 0.001 from right to left, respectively. The left panel is for a galaxy survey with / < 25.6 
and 0.4 < z < 4.2 and the right panel is for a "toy" galaxy survey capable of observing 10 -3 galaxies per (Mpc/h) 3 in 
the range 4.2 < z < 8.0. Comparing these plots we see that the former survey does a better job of constraining r but 
the difference between these surveys diminishes for smaller values of A p . We conclude that surveys that can observe 
higher density of galaxies are preferable to the ones that go to the larger redshifts but observe relatively low galaxy 
surface density. We can also see from these plots that the contours of constant ay become steep at small A p so that 
a small improvement in the detectors sensitivity goes a longer way towards detecting GWs than does decreasing the 
beam width. 
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FIG. 3: The error bars for measuring the tensor-to-scalar ratio r by constructing either B l LM (red circles) or E l LM (green 
crosses). The galaxy survey has the specification J max = 25.6 and 0.4 < z < 4.2. The B polarization is more powerful because 
it is not contaminated by scalar modes. 
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FIG. 4: To find out how the signal for detecting r is distributed in multipole space we have plotted the error on measuring the 
tensor-to-scale ratio versus the maximum I2 in the sum in Eq. [30] . We see that up to the i™ 3 * ~ 4000 the signal increases with 
increasing Zf" 1 * after which it levels off. 



IV. SYSTEMATIC EFFECTS 



An understanding of systematic effects is crucial to properly interpret any data set. Biases can be introduced to 
our parameter estimations either by unaccounted for instrumental errors or by contaminating signals, such as residual 
foregrounds. One way of finding the extent to which the data has been contaminated by systematic effects is by 
identifying general properties that the data must satisfy, e.g. invariance under parity, where breakdown of any of 
these properties is considered a strong indication of an spurious signal. Another way of quantifying the biases is by 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

A,, (^/K-arcmin) A p (//K-arcmin) 



FIG. 5: Contour plot of the estimated error in the tensor-to-scalar ratio, oy, as a function of detector noise and beam width of 
a CMB probe. The contours from right to left are for oy = 0.1, 0.05, 0.01 and 0.001, respectively. The left panel is for a galaxy 
survey with maximum limiting magnitude in the / band of J ma x = 25.6 and redshift range of 0.4 < z < 4.2. The right panel 
is for a survey capable of observing a constant density of 10 -3 galaxies per (Mpc/h) 3 from redshift z = 0.4 up to z = 8. The 
latter survey has a smaller number density at low redshifts (z < 3.5) than the former survey but it has a non-negligible surface 
density of ~ 10 6 galaxies per square radian up to the highest redshifts. 



calculating their effects on parameter estimations, either analytically or by running a suite of Monte-Carlo simulations. 

In this section we estimate biases from the gravitational lensing of the CMB, polarized point sources, and incomplete 
sky coverage. These effects bias our estimators because all of them can introduce correlations between the observed 
CMB polarization and large scale structure for non-equal values of I. This correlation is exactly what we have exploited 
in our method to estimate the GW signal, so the extra correlation from lensing and point sources, if unaccounted 
for, can be incorrectly interpreted as a GW signal. It will be shown that these foregrounds will not add any bias 
to our estimator for r and only increase the error bars by adding to the cosmic variance. In the case of incomplete 
sky coverage we notice that it is no longer possible to work in multipole space as the E and B modes are non-local 
functions of the polarization on the sky and therefore cannot be calculated unambiguously. Given that most of our 
signal comes from very large angular scales (L = 2) the cuts will have a dramatic effect on our estimates of r. Here 
we set out to calculate these effects quantitatively. 



A. Weak lensing of the CMB 

As mentioned before, the observed CMB polarization, X ; °^ s , where X can be either B or E, can be written as a 
sum of the polarization produced around the time of recombination, Xf^° m , during the reionization era, -XJ™ * 1 , from 
the gravitational lensing of the primordial polarization, Xj^ s , from the polarized point sources, Xf^' , from polarized 
galactic emission, Xf^ , and finally from detector noise Xf™ c : 

vrccom i v~rcion . -^lcns , -v^p.s. . -v^gal . xrnoise f ACi\ 

X lm - A im + A Zm + A im + A *m + A im + A im ■ l 4U i 



The polarization produced by the gravitational lensing of the primordial polarization can be written as 22[ : 



Ellens \ ^ \ ^ j / -] \mi / ^1 ^6 ^5 \ T^lcns t^i 

E im = 2^ 1^ W6m.(-1) I mi _ TOR _ TOr ) , 2F hhh a hh i a E l 

h h I 



rccom 



B lm S = ^n 5 (-l) mi ( ^ __^ g ) 2F/ 1 ?/i 6 7hM6K°m 6 m > ( 41 ) 
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where we have ignored the primordial B mode polarization. Here 2-F/ 1 e ; Ils ; is defined as 



- Hh + 1) + + 1) - h(h + 1)] y ' 2 '^ 1 " 2 '^ 1 " 2 '^" ( J J 1« 2 ) (42) 
and <fr is the lensing potential, which can be written as a weighted line of sight integral of the gravitational potential, 

* (i)= - 2 r %°t4 * |a ' wiA <43) 

where tjq is the conformal time now and rji s that at the last scattering surface. 

The bias added to our estimator from the weak lensing of the CMB can then be found by calculating (j3 l LM ^J. A 
straightforward calculation gives: 

Bl M (len S )) = ^EEt-W^ ™ 2 _M)^f 2 (i)(7W 2i fcAU^ + aw 2 LteAL m2 >) 

l\m\ l 2 m,2 

= 0. (44) 



ing to the background covariance matrix (B™ M Bl, M , ), i.e. it adds to the cosmic variance. To calculate this noise, 



Therefore, the lensing does not add any biases to our estimator. However, lensing acts as a source of noise by contribut- 

^LM D L'M' 

we need to find ensemble averages of the form (ij* A[* m! 1'°^ Aj 4mi ) . Ignoring non-Gaussianities, this ensemble 
average can be expanded into three terms. The noise coming from the terms of the form (A'° i ^*X'° 3 b ^ l3 )(A;* m2 A| 4m4 ) 
have already been taken into account by using the total (j EE '° hs a nd (j BB > oh * j including lensing, in our formulas for 
the noise covariance matrix (see Eqs. I16||2T)1) . The extra noise coming from the other two contractions is: 

(Bt M Bi, M ,}(lens - lens) - 6 W 6 MM , HJ?L 

EEEEK- 1 )' 1 " 13 ^ ^ 5 2 & < > { t \\ i}c^ ccom c^cf;\ (45) 

h h h U h L J 

where I > is the Wigner 6-J symbol, and < I2I4J5L >= 1 when I2 + h + h + L =odd and zero otherwise. The 



correlation between the lensing potential and the galaxy distribution at redshift slice i is defined as (4>i m A\, m ,) = 
Cf l Su> 8 mm ' ■ However, this term is subdominant compared to the term we have already considered. 

To find out to what extent weak lensing of the CMB contaminates the signal we have computed the expected errors 
on r assuming no lensing of the CMB. Even though this is an unrealistic assumption it will serve to illustrate the 
degradation in the signal caused by the gravitational lensing of the CMB. We find that the error on r reduces to 
0> ~ 0.04 for complete cleaning of the lensing, compared to ay ~ 0.09 for no lensing cleaning, so at best the signal 
can be improved by a factor of ~ 2. 



B. Polarized Point Sources 



Since the distribution of polarized point sources follows that of dark matter they can add spurious signal to our 
estimate of r, if uncleaned, by correlating the galaxy distribution and the CMB polarization fluctuations. In this 
section we quantify this effect by calculating the bias and noise introduced to our estimator of the tensor-to-scalar 
ratio by polarized point sources . 

The E/B polarization multipoles for point sources can be found to be: 

Efm±iBl^ = J d 2 h ±2 Y l * m (n)[Q^-(n)±iU^-(n)} 
- J d 2 n ±2 ^ m (n) 
= E ±2^ m ( ni )5 ie ±2 ^. (46) 



. i 
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Here Si is the polarized flux of the ith point source, hi is its direction on the sky, and is the direction of its 
polarization with respect to the eg unit vector on the sphere. 

Similar to the lensing case, to find the bias introduced to our estimator from point sources we need to calculate the 
contribution to our estimator (b lm ^J from point sources, where now the ensemble average also includes averaging 

over Si, hi, and uii. This average is trivially zero because {e ±2uJi ) [Aji = 0. So the polarized point sources do not add 
any bias to our estimator. 

Calculation of the noise contributed by point sources is a bit more involved. The final result is: 



A*' B A j ' B (T S 2 ) 

{B? M Bl, M ,){v. S . - p.s.) = 5 LL ,5 MM , ^± ^ p} 
E E E E E(-D' 3+ ' 4+i5 ^ (L)9^(Lh hlML B^ u 



h h h 
-2 2 



h h L 
h h h 



(47) 



where the bispectrum B is defined as: 



\ (2TO2 ferns I. 



h h U 
vri2 ms 7714. 



Here A p s - is the projected overdensity of point sources on the sky. 



B 



/ _>/ ". 1 1 



(48) 



C. Incomplete Sky Coverage 

A real CMBxLSS cross-correlation will have only partial sky coverage, due at the very least to the Galactic Plane, 
and possibly to additional observing constraints (e.g. hemispherical). This is troublesome for our method because 
most of the signal comes from large angular scales (see Fig. [3]). Even though it is possible to clean these foregrounds 
by, for example, using the frequency dependence of synchrotron and dust emission, this procedure might introduce 
unwanted correlations in the residual polarization maps, which in turn can be mistaken for a signal in our estimator. 
In this paper we take a more approach whereby we completely cut a portion of sky that is most severely contaminated 
by polarized galactic emission. However, by doing this we can no longer use our estimators, written in multipole 
space, since the Ei m and Bi rn coefficients are nonlocal functions of the polarization in real space and for a given I 

they cannot be reliably estimated for all possible m. Consequently, our estimators for Elm Blm are not applicable 
in this situation. 

Here we outline the procedure to overcome this issue and refer the interested reader to Appendix [C] for details. To 
solve this problem, one needs to work in configuration space. We first find an estimator for the average polarization 
generated at a given redshift bin "i" coming from a given direction n on the sky, X l (n), where \ stands for the Stokes 

parameters Q or U. We then compute the 2N z N p - lx x 2N z N p - lx covariance matrix C^ AXB = ( n ^)Xi3 ( n s)^ for 

data consisting of Q l (n) and £/*(n) measured in iVpi x pixels and N z redshift bins. This covariance matrix consists 
of a part proportional to r coming from gravitational waves and a part from all other sources. In a Fisher analysis 
these modes can be projected out by giving them formally infinite power. Finally, we calculate the error bars on r by 
computing the Fisher matrix using a Monte Carlo code. 

We report our results in terms of a "degradation factor" /dcg defined as the ratio of oy for a cut sky to its value for 
complete sky coverage. We have computed /dog using the Monte Carlo trace code of Ref. (23|. Computations were 
done pixelized at HEALPix resolution 8, which has 786 432 pixels |24[, although convergence was found at resolution 
7. We chose for the E-mode projection a Gaussian prior X( — 50e~ e ^ e+1 ^ a I" 1 with a = 0.079 rad = 4.5 deg; the 
important point is for the low is to have 3> 1. A total of 200 Monte Carlo trace realizations were run for each sky 
cut. We find a degradation factor of /d cg = 0.32 ±0.01 for a sky cut of |6| > 10° (/ sky = 0.83) and / dcg = 0.056 ±0.004 
for a sky cut of \b\ > 20° (/ s k y = 0.66). Note in particular that /d og < /sky for these cases, and that /d cg ~ /sky is a 
poor approximation. 



V. DISCUSSION AND CONCLUSION 



Part of the polarization pattern of the CMB is generated by the Thompson scattering of photons off of free electrons. 
The strength of this effect is proportional to the local quadrupole moment of the radiation and the number density of 
free electrons. The polarization that is generated during the reionization era is therefore modulated on small angular 
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scales by the varying electron number density. This variation in electron number density itself is correlated with the 
distribution of galaxies. So, by looking for correlations between the small scale CMB polarization fluctuations with 
the galaxy number density at a given redshift one can determine the local quadrupole moment of the CMB at that 
redshift. These quadrupoles at different patches of the sky and at different redshifts then provide us with a map of 
the quadrupole moments during the reionization era .39]. A small part of this quadrupole pattern can be produced 
by the tensor modes of fluctuations, i.e. gravitational waves. Therefore, the correlation between galaxy distribution 
and the CMB polarization anisotropies can be used to constrain the strength of the primordial GWs, which are of 
the utmost importance for physicists. They provide us with information about the physical processes at work in the 
early Universe with energies way beyond the reach of any terrestrial experiment. 

In this paper we have investigated the prospect of using this method to measure the tensor-to-scalar ratio r, which 
is a measure of the strength of the GWs. We have constructed a full sky estimator for measuring the average B 
and E polarization generated at any given redshift bin during the reionization era by using the observed distribution 
of galaxies at that redshift and the observed CMB polarization pattern. The result of this exercise was then used, 
together with linearized perturbation theory and Fisher formalism, to predict the prospects of detecting a GW signal 
from future experiments. We have found that a CMB experiment with noise parameters of Ap = 2/iK-arcmin and 
0fwhm = 2 arcmin, together with a LSST-like galaxy survey [40] | with a limiting magnitude of I < 25.6 and a similar 
northern component 4l[ can be used to constrain r to oy ~ 0.09. Even though these are photometric redshift surveys 
the small fraction of incorrect redshift determination will not be a source of problem for our method because it cannot 
produce spurious B modes in P. 

To reach this level of accuracy many obstacles need to be overcome, including the proper handling of the systematic 
effects. Here, we investigated two of these effects: the polarization induced by the weak gravitational lensing of 
the CMB and polarized point sources. We found that neither of these contaminants adds a bias to our estimator. 
However, they increase the error bars on r by increasing the cosmic variance by adding to the observed CMB power 
fjBB,ohs an( j (jEE,ohs ^ showed that lensing increases ay by a factor of ~ 2, so lensing cleaning techniques might 
partially help to improve the signal. The most important factor, however, turns out to be the incomplete sky coverage 
where the degradation in signal to noise is more severe than the naive estimate oc /sky For example, for a sky 

cut of \b\ > 10° (/ sky = 0.83)) the depredation factor is / deg = 0.32 and for \b\ > 20° (/ sky = 0.66) it is / dcg = 0.056. 



Acknowledgments 

E.A. and C.H. were supported by the U.S. Department of Energy (DE-FG03-92-ER40701). C.H. is also supported 
by the National Science Foundation (AST-0807337) and David & Lucile Packard Foundation. E.A. acknowledges 
support from the National Science Foundation (grant AST 07-08849) during part of this work. He thanks Laura Book 
and Ben Wandelt for their help and Avi Loeb and John Kovac for useful conversations. 



Appendix A: Derivation of the quadratic estimator 

In this section we derive formula [14] for the E l estimator. The derivation for B l is similar. This calculation closely 
resembles the lensing reconstruction method of Ref. [l2[ • 
We start from a general quadratic estimator of the form 

EIm ^fEE iit-M) + 9i'l B (L)BfX) (Al) 



and our goal is to find the unknowns Al, g EE and g EB . To find them we use two common criteria for a good 
estimator, that is that it should be unbiased and have minimum variance. 
For the estimator to be unbiased we need 

(EIm) „ , =Ei M , (A2) 

small angles 

that is the ensemble average of the estimator over different realizations of the small angular scales of the polarization 
and galaxy distribution should be equal to its true value. We must then have 

,E ( 1 \M I h h L \ ( ' J,EE i T \i pobs \i,obs\ , J,EB / T \/ D obs 



13 



The correlations between the observed galaxy and polarization distributions can be calculated using Eqs. (122II23[) 



l\mi 



l h h 

—m mi mi 

I h V 

— m m\ —to' 



(J9 l 9 % 

Flilh{<XlihlElimi - lhhl B hrni) Ti 5i 2 ,i>8 m2 -m' (-l) r 



/ F hW (a h w E l hmi - 7i 1 M'-B; imi ) 



b> 



(A4) 



Using the identity 



h h h 
-mi — m 2 — TO3 



the fact that an odd permutation of columns 



) V 7711 7712 m3 

of the 3-j symbol introduces a similar factor of (— l)'i+'2+Js ) the fact that for the 3-j symbol to be non-zero we must 
have mi + mi + 7713 = 0, and changing the indices from (1%, mi) to (L, M), we finally find 



obs a i,obs\ 



small angles 



M 



LM 



I I' L 
m to' — M 



fill' ( a Lii'E l LM - JLU'B 1 LM ) , 



flw = Fuv ■ 



cf, 



A similar calculation gives 



/ nobs a i,obs\ 



small angles 



Em: 



A/ 



LM 



? V L 

TO to' —A/ 



flw \lLii'E l LM + a LW B LM ) 



(A5) 



(A6) 



Plugging these results into Eq. (|A3|) we find a relation between the unknown coefficients of the form 

9hi 2 \ h ) _ 9i x i 2 \ L ) _ i,E tT ^ 

— = 9r,iA^)j 

a Ul i 2 luih 

A B = (2L+l)feft hh gffi(L)\ . 



(A7) 



minimum variance 



Using these relations, the only unknown coefficients are g\^ 2 {L). To find them we demand that the estimator have 

(A8) 



d ( E L*M E LI 



J LM LM I 



= 0, 



where the variance can be calculated to be 



K^LM LMl — °L 



E l E l 



{A 



i.E\2 



2L 



^( 9 ;|(L)) 2 < ll2 c^'» bs , 



hl 2 L 



(\ 1 2 ^ EE, obs . 1 a r^BB, o\n 

[\aLhl 2 \ C h +\1L1 1 1 2 \ C h ■ 



Using this in Eq. (|A8|) we finally obtain 



flh 



E syg'g*, obe' 



^^l\liL l 2 



(A9) 
(A10) 

(All) 



Appendix B: Signal covariance matrix 

In this section we calculate the signal covariance matrix Cf B . We start from the line of sight solution to the 
Boltzmann equation for i?-type polarization, see e.g. (25|, 

E\ m (k) = E\ m { V = m , k) = -VE(2l + 1) / d V g(v)Pm(V, k)e lm (k( m - V )). (Bl) 
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In coordinates where z || k, m = 0, 1, 2 correspond to the scalar, vector, and tensor modes, respectively. The projection 
functions for scalar and tensor modes are 



e i,o( x ) 



I3 (l + 2)lj l (x) 
8 (Z — 2)! x 2 ' 

-Ji W + J; W + —jT~ + — — 



and the source functions are 



p — 

1 m — 



1 

To 



e 2 



VEe 2 



10 



-e 



2,m ; 



(B2) 



(B3) 



where we have ignored the polarization term in the source equation above compared to the temperature term. The 
temperature multipoles themselves can also be found as line of sight integrals 



Q;,o(?7, k) 
21 + 1 

2Z + 1 



dr\ e 



' P -r(n,n') 



(t-00,0 + r\t + * - $)jj,o,o + rvoji,ifl + fP ji i2 ,o 



dj/e^ [tP 2 (t/, fc) - ij(+ 2 )(r/, ji, 2 , 2 (k( V - if), 



(B4) 



for the scalar and tensor modes, respectively. Here the optical depth is r(r], r/) = C f(rj")dr]" and H^ +2 ^ is the 
amplitude of the right handed gravitational wave. The projection functions are 



ji,o,o( x ) = 3l( x )> 3i,i,o = fi( x ), 3i,2,o( x ) = r [3j"(x) + ji(x)] , 



JJ,2,2(ar) = 



'3 (l + 2)lj,(aQ 
8(1-2)! x 2 



(B5) 



Assuming that recombination occurred instantaneously at r\ = 77* , neglecting the scattering source term rP 2 compared 
to redshifting H, and writing the metric perturbation in terms of the intial metric perturbation and the transfer func- 
tion H^ +2 \rj.k) = T(rjk)H^ 2 \k), we find for the temperature quadrupole coming from right-handed gravitational 
waves 



e 2 , 2 (r7, k) = -5H { + 2 \k) [I(k V ) - J(fofc)] 
I{x) = 



dx . j ]2,2,2{X - X ). 



dx' 



The transfer function in the matter dominated universe is 



T{x) 



3ji(z) 



(B6) 
(B7) 

(B8) 



We show a plot of I(x) in Fig. [5] To find Cf E we notice that the multipoles Ej (k) in equation IB 1 1 are different 
from the Ej m (k) we used in our calculation in the main text because instead of ± 2 Y l TO (n) they are expanded in terms 
of 



(-0' 



47T 

27+1 



±2^7 m (n). 



(B9) 



Correcting for this normalization and assuming equal amounts of left and right circularly polarized waves, the power 
spectrum of the tensor modes can be found as 



n E i E j 



7T(2/ + l) 2 



dkk'P, 



2 r>E l E ] 



1,2 



(AO, 



(BIO) 



where P E 2 E (k) is defined as 



(Er 2 (k)EUk')} = {2^fpg E \k)5 3 {k k') 



(Bll) 
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FIG. 6: [Left] The evolution of the temperature quadrupole generated by the gravitational waves. [Right] Transfer function 
S l £ vs wavenumber k, shown for different choices of multipoles L and of redshift bins "i". 

Using equations EU |H EU EQ] and [BIT] one finally finds 



C, 



dfcfc 2 5;" B (fc)s*/" E (fc)p H ( +2 , (fc), 



= ~Y / '%! '/>'/. j!/.'!//,, -- //)) (./(/■•//) ' /!/■//; !! . 



(B12) 



where P H (+$) is the power spectrum of the initial right handed GWs. To make a connection to the definition of Ph in 



reference (T(| we notice that H^ 2 "* = —\\ \{h+ — ih x ) and therefore 



1 „r, 



P H ^(k) = -((\h + \ 2 ) + (\h^)) = -(\h\') = 



Ph(k) 



(B13) 



An equation similar to Eg. IB12l can be found for Cf T B] by substituting ei t % — > /3i % where 



(B14) 



In the right panel of Fig. |6] we show the transfer function S(k) for the quadrupole moment of the B polarization 
generated at the first redshift bin. 



Appendix C: Sky Cuts 



Here we present the detailed calculation the degradation in our estimation for the tensor-to-scalar ratio caused by 
the sky cuts. As mentioned before, these cuts are troublesome for a multipole-space based estimator because E and 
B modes cannot be found unambiguously in this case. Fortunately, the reconstruction of the mean polarization field 
in the i th redshift slice, P l (n), is a "local" operation in the sense that it only depends on information within a few 
degrees of n (see below and Ref. [Hj]). However, this still leaves us with the problem of computing the large-scale 
power spectrum C e lJ and forecasting its errors in the presence of a cut sky. This is the same problem that occurs 
in the context of the CMB power spectrum and many works have been devoted to the problem [2a - l29| . 

If we define the filtered fields 



-P F (n) 




±2^ ljmi (n), 



Z 2 m 2 W *°/ 2 



(~<9 a 

l 2 A»,obs y (~\ 

g i g i ,obs hm 2 x (2,m 2 V i V! 
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then one can prove that 

± P F -(n) = ± P*(A)A*'<(A) = £ ElM ^.f LM ±2Y LM (n). (C2) 

IM L 

We have then accomplished our goal of constructing an estimator for the average polarization generated at a given 
redshift bin that is not affected by the eliminated portions of the sky, since the quantities ±P F (h) and A F,l (n) are 
local in configuration space. We can further simplify the analysis by noting that the A l £ B are very nearly independent 
of scale for small L, where most of the signal comes from. Thus we can pull this factor out of the sum and find an 
estimator for the non-filtered average polarization generated at the redshift bin "i" 

± P\h) = £ (El M ± iBl M ) ± 2 Y L>M (n) = A l f ± P F <\h). (C3) 

LAI 

The estimators for the Stokes parameters, Q l (n) and f/ l (n), can then be found as 

Q\n) = + P ' (n) + - Pt(n) = Y / E i LM Y L E M (n,Q)+S i LM Y B M (h,Q), 

LM 

U*(A) = +Pt(n) 2l ' F ' (n) = E ^lmYlm^, U) + B^Y^in, U), (C4) 

1 LM 

where we have defined 

Y F M (h,Q) - Y B M {n,U) = \( 2 Y LM {n) + - 2 Y L<M (h)) , 

Y B M (h, Q) = -Ff M (n, U)= l -{ 2 Y LM (n) - -2^(6)) • (C5) 

We write down a vector x consisting of the mean polarization Q'(n^) and U % (nA) in the A th pixel, of length 2N z N p i x . 
Then x has a 2N z N p i x x 2N z N p i x covariance matrix 

C% AXB (n A ,n B ), (C6) 

where xa indicates a choice of Stokes parameter (Q or U). This contains a contribution Co arising from all sources 
other than gravitational waves (including the scalar post-reionization scattering signal and noise), as well as a contri- 
bution associated with tensor modes, 

C = C + rC r . (C7) 

If it is desired to remove certain modes from the data (e.g. E- modes), then they can be incorporated in Co by 
introducing such modes with formally infinite power. (This is a common foreground template projection method in 
CMB data analysis, e.g. [30}.) 

Under the null hypothesis of no primordial gravitational waves, and assuming Gaussian signal and noise, one can 
write down the Fisher information, 

F rr = -Tr (C 1 C r C 1 C r ) , (C8) 

— 1/2 

and forecast an uncertainty oy = F rr ' . The assumption of Gaussian signal and noise is expected to be valid in this 
case: the signal (primordial gravitational waves) are Gaussian and feed linearly into P*(fi). The noise (associated 
with the product of primary CMB and large scale structure) is not Gaussian, but its average in a large patch of sky 
(tens of degrees) contains contributions from many arcminute-scale patches, so we expect the noise contribution to 
P l (n) to be Gaussianized by the central limit theorem. 

The actual computation of F rr is simplified if we consider only the i?-mode quadrupole, 1 = 2, and note that the 
noise power spectrum is white at low I (i.e. it becomes independent of I and equal for E and B modes). In this case, 
the matrices have the form 

2 

[C r ]^ AXB (h A ,h B ) = Cf BJ £ Yi:(n A , XA )Y 2 B m (h B>XB ), (C9) 

m=-2 
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where Y 2 B n (n, x) denotes the X component of the B-mode tensor spherical harmonic Q evaluated at n (see Eg. lC5|) . If 
we insist on projecting out all -E-mode signals (so that there is no possible contamination from the scalar reionization 
signal), then the zero-tensor covariance matrix is 

[C ]^ AXB (n A ,n B ) = N^5 XAXB d^\n A - n B ) 

+ E A ^ E Y 2 E n :(n A , XA )Y 2 E m {h B , XB ), (CIO) 

1=2 m=-l 

where A™ is the white noise level (still a matrix because it depends on the redshift slices), and the second term 
projects out i?-mode signals. The projection parameters A^ should formally be taken to oo. 

Under the above assumptions, the matrix inversions and trace factor into pieces that depend on the sky coverage, 
and pieces that depend on the signal and noise power spectra (C B B ° and A™): 



F rr = ^Tr (n -1 C| b N- 1 C§ Tr (K _1 SK _1 S) , (Cll) 



where the first trace is of a product of N z x N z matrices, and the second trace is of a product of the 2iV pix x 27V p j x 
matrices: 

K XAXB (n Al h B ) = S XAXB S^(n A - n B ) + g A, £ Y 2 E *{h A , XA )Y E m {h B , X b) (C12) 

1=2 m=-£ 

(again, valid only in the Xe — >• oo limit) and 

2 

S X axb (»a, n B ) = E Y 2m{n A , XA )Y 2 B m (iL B ,XB)- (C13) 

m = -2 

In Eq. ([Clip , the first trace depends on detailed galaxy and CMB properties and has been computed in the main 
paper. The second trace enables us to define a "degradation factor:" 

Tr (K^SK^S) (true) 
/dog = T7(K rT SK rT S) (all sky)' (C14) 

— 1/2 

It is jdeg that this appendix aims to compute. For a cut sky, the uncertainty on r increases by a factor of / deg . The 
analytic expectation is that in the harmonic basis and on the full sky, K becomes a diagonal matrix with Is in the 
B-mode positions and oos in the iS-mode positions (because of the formally infinite iJ-mode terms added); S should 
have Is in the diagonal 1 = 1 £?-mode positions and Os elsewhere. Thus the denominator should simply be the number 
of £ — 2 B-modes, i.e. 5. We find that numerically this is indeed the case. 
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